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By resorting to recent results on the relativistic currents for mixed (flavor) fields, we calculate 
a space-time dependent neutrino oscillation formula in Quantum Field Theory. Our formulation 
\ provides an alternative to existing approaches for the derivation of space dependent oscillation 

. formulas and it also accounts for the corrections due to the non-trivial nature of the flavor vacuum. 

' By exploring different limits of our formula, we recover already known results. We study in detail 

^S) ' the case of one- dimensional propagation with gaussian wave-packets both in the relativistic and in 

the non-relativistic regions: in the last case, numerical evaluations of our result show significant 
deviations from the standard formula. 
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■ I. INTRODUCTION 

> " 

, The recent experimental evidence of neutrino mixing and oscillations is the first clear sign of physics 

' beyond the Standard Model 7] . Thus much effort is currently devoted to the full understanding of such phenomenon, 
\ from the issue of its origin to more phenomenological ones. 

' In this framework, there has been recently a remarkable prog ress in the study of the problem of field mixing in 
^ ; Quantum Field Theory (QFT) [ili|li|ll|lllillil|p,|3lHlIllIllp-Iti^ important to remark that, in spite 
. of the fact of being usually treated in Quantum Mechanics (QM), the mixing between states of different masses is not 
' even allowed in a non-relativistic theory due to the Bargmann su per selection rules . The problem of constructing 
O-i, an Hilbert space for flavor states is indeed a long standing one |22j which has even been claimed to be impossible 
O •■ to achieve ^2§\ (see however Ref.p3| for a criticism of that argument). As first shown in Ref.l3|, the difficulty lies in 
1^ , the fact that the Hilbert spaces for particles with definite flavor and those with definite mass are actually orthogonal 
• in the infinite volume limit. This result has been subsequently generalized for any number of generations and for 
^ [ different types (fermionic or bosonic) of fields [lol IT^ IT^ . It has also emerged that the use of the correct Hilbert space 
(i.e. the one for the flavor fields) leads to corrections to the usual Pontecorvo oscillation formula. 

The necessity for a full QFT treatment of mixing also stems from a more phenomenological perspective, namely 
from the calculation of (flavor) oscillation formulas. The usual treatment gives indeed in a very simple way an 
oscillation formula in time, which is actually not of great use when discussing current experiments where the distance 
source-detector is measured rather than time (data being collected over large time intervals) 24]. In order to derive an 
oscillation formula in space, one usually converts the time oscillation formula by means of some assumptions (equal time 
assumption, classical propagation, etc.) which are however questionable in many respects (see Ref. 25] for a review). 
Thus several QFT approaches have been developed |2^ I2M I29I [sS] , such as the external wavepacket models or the 
stationary boundary conditions models, which attempt a more realistic description of neutrino oscillations. It seems 
indeed that the main features of such formula are now quite well understood 25], including the concepts of coherence 
length, localization terms, dispersion times, etc. 

On the other hand, all these calculation schemes have been developed for highly relativistic neutrinos and/or for 
nearly degenerate masses and in most cases the spin structure of the wavefunctions was neglected (i.e. neutrinos were 
treated as scalars). Also, the identification of the Hilbert spaces for mass and flavor neutrinos was implicitly assumed 
in the choice of the propagators (see Ref.Q for a discussion of flavor propagators). 

A completely different approach to the problem is indeed possible and was advocated in Refs. fsills^ : an oscillation 
formula in space could be derived in a straightforward way, by considering the flux of neutrinos through the detector 
as the integral over the measurement time and detector surface of the relativistic flavor current for the oscillating 
neutrinos. An attempt in this direction is contained in Ref. 'Sl], which however deals with QM and therefore uses a 
probability current derived from Schrodinger equation. As pointed out in Ref. [32j] . the obstacle to the extension of 
this result to QFT lies in the difficulty of constructing a relativistic current for mixed particles, which boils down 
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again to the problem of the definition of the Hilbert space for such (not on-shell) particles. 

The solution to this problem was recently given in Ref.|^3- it turned out that it is indeed possible to define these 
(non-conserved) flavor currents in a consistent way by use of previous results on the flavor Hilbert space. An analysis 
of the currents associated to mixing for the case of three flavor neutrino mixing with CP violation and for boson 
mixing can be found in Refs.^lft. .llj. 

In this paper, we calculate the oscillation formula by using the relativistic flavor currents above mentioned. For 
simplicity, we only treat the case of mixing among two generations of Dirac neutrinos: the extension to Majorana 
fields and to three flavors will be given elsewhere. The formula we obtain in a very direct way has the full space-time 
dependence and contains previous results obtained by use of the flavor charges. We then obtain a general expression 
for the electron neutrino flux in three dimensions with spherical symmetry. Although a fully three dimensional 
analysis is possible within our formulation, we study in detail the one-dimensional case with gaussian wavcpackets, 
which is indeed the one most frequently considered in literature. We show how in the highly relativistic limit and by 
integrating over an infinite time, the standard space dependent oscillation formula |25l 1261] is obtained, which includes 
the coherence length and the localization term. However our formula is also valid for non-relativistic neutrinos and 
it accounts for flavor vacuum effects. Numerical evaluations of our result show deviations from the usual formula in 
the non-relativistic regime. 

The paper is organized as follows: in Section H we present a derivation of the flavor currents in the line of Ref.|l3| 
and give some details on the construction of the flavor Hilbert space. In Section III we first obtain the expectation 
value of the fiavor current in the most general case and then a more explicit expression in the case of a spherically 
symmetric emission. Then, in Section IV, we specialize to the one-dimensional case and obtain the space dependent 
expression for the electron neutrino flux. Further analysis is done by choosing gaussian wavepackets. Section V is 
devoted to conclusions. 



II. RELATIVISTIC CURRENTS FOR MIXED FIELDS 

Following Ref.|l3j|. let us consider the following Lagrangian density describing two Dirac fields with a mixed mass 
term: 

= ^fix){t^-M)^f{x), (1) 

The mixing transformations 

/ cos 6 sinO \ , 
- [sine cos0j*™(")' (2) 

with being the mixing angle and 'I'^ = {1^1,1^2), diagonalize the quadratic form of Eq.(^ to the Lagrangian for two 
free Dirac fields, with masses mi and m2'. 

£{x) = ^ - Md) *™(x) , (3) 

where = diag{mi, m2). One also has rtie — mi cos^ 9 + TO2 sin^ 9 , = mi sin^ 9 + TO2 cos^ 9 , = (m2 — 
TOi) sin 6* cos 6* . We assume m2 > mi and 9 in [0,|^]. 

The Lagrangian Eq.(^ is invariant under global U{\) phase transformations of the type = e'" 5'^: as a result, 
we have the conservation of the Noether charge Q — J d^xl'^{x) (with /^(x) = ^m{x)'y^ ^'^(a;)) which is indeed the 
total charge of the system (i.e. the total lepton number). 

Consider now the SU{2) transformations acting on "^m- 

^'^{x) = e^'^^^^ ^^{x) , J = 1,2,3, (4) 

with aj real constants, Tj = crj/2 and aj being the Pauli matrices. 

For mi m2, the Lagrangian is not generally invariant under Q and we obtain, by use of the equations of motion, 

SC{x) = *m(x) [r,-, Md] ^'m(x) = -^ct] d^J^ .j{x), 

^™,,(-^) = ^mix) r r, *„.(a;) , j - 1, 2, 3. (5) 

The above analysis is valid at classical level. We now quantize the free fields vi and V2 in the usual way (see 
Appendix A for conventions). 



C{x) 

where — {i^e, v^) and M = 



me mefj. 

m, 



■^f{x) 
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Then the charge operators Qmjit) = J d^x J^j(a;), satisfy the su{2) algebra [Qmj(t), Qm,fc W] = i ejki Qm,ii^)- 
The Casimir operator is proportional to the total charge. From © we also see that Qm,3 is conserved as M^i is 
diagonal. Let us define the combinations: 

Ql = + Qrn,^ , Q2 = - Qm,3 (6) 

Q^ = Y. f ^'1^ - , * = 1,2, (7) 

r •' 

where the last expression has been normal ordered, as usual. These are nothing but the Noether charges associated 
with the non- interacting fields vi and V2'. in the absence of mixing, they are the flavor charges, separately conserved 
for each generation. Observe now that the transformation 

is just the mixing Eq.Q. Thus Geit) = e^^^Q^.^C*) is the generator of the mixing transformations (see Appendix 
A). In Ref.jQi], it has been shown that the action of Ge{t) on the vacuum state |0)i.2 results in a new vector (flavor 
vacuum) |0(i))e,u^ G^^(i)|0)i_2, orthogonal to |0)i,2 in the infinite volume limit. 

Following Ref. , in accordance with Eq. © , we define the flavor charges for mixed fields as 

Q„{t) EE Gg\t)Q,Ge{t) (9) 

with ((T, i) ~ (e, 1), (/i, 2) and Qe{t) + Qi_i{t) = Q. They have a simple expression in terms of the flavor ladder operators: 

QAt) = / d'^4{x)Mx) = E/'^'k {all{t)alAt) - /3l,.(i)/3:k,.W) , ^ = e,M. (10) 

These charge operators then act on the flavor Hilbert space Ti-cf^ (see Appendix A). For the single neutrino and 
antineutrino states of definite momentum and helicity^: 

l^.(t)) ^ a;,|,(t) |0(t))e,^ ; M))^ PlUt)m))e,,, (11) 

one obtains 

QAt)K{t)) = WAt)) ; Qa{t)\9At)) = ~\Mt)) ■ (12) 

In the following, we will use the flavor currents: 

J^{x) EE D^{x)^'',y,{x) = G-\t)J>^{x)Gg{t), (13) 
with (cr, i) ~ (e, 1), (/i, 2). From the global U{1) invariance follows the continuity equation: 

d^[Jlt{x) + Ji^{x)\ = 0. (14) 



III. SPACE-TIME DEPENDENT NEUTRINO OSCILLATION FORMULA 

In this Section, we consider the calculation of the general (i.e. space-time dependent) neutrino oscillation formula, 
by use of the flavor currents introduced in §11. In Section IV, we then extract from this the space dependent formula 
by integrating over time. 

As already discussed in Ref.[3ll Is^]. although in the context of QM, what one is ultimately interested in current 
experiments is actually the flux of neutrinos of a given flavor through the surface of the detector in a (large) 
measurement time T. In the case of electron neutrinos, this quantity is given by: 

$,^^,^(L)= / dt ( {u,\Jl{^,t)\u,) dS\ (15) 

J to Jo. 

where L is the distance source-detector. Our aim is now to calculate the expectation value of the flavor four-current 
density on a localized neutrino state, defined as a wave-packet in flavor Hilbert space. 



^ Note that there is no ambiguity here in defining the flavor neutrino states with definite momentum. 
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A. Expectation value of the flavor current density 



We define an initial^, localized state |i/e(xo,to)) with definite flavor e : 

where ^ is a normalization constant, /(k) is the form of the wave-packet and |0)e_^ the flavor vacuum at time t — tg. 
This corresponds to an electron neutrino being emitted at (tg, Xq). For convenience we set (to, Xo) = (0, 0, 0, 0) in the 
following. Note that the wave packet is deflned in the Hilbert space of flavor fields. 
The normalization of this state is given by: 



1 = (lyeWe) = \A\' J d'k J d'pr{k)f{p){alM,c^l]em = \A\' / d^k |/(k)| 
From Section II, the electron neutrino four-current is given as 



(17) 



where F'^ = 7^7'^. In the following we will use the chiral representation of the Dirac matrices (see Appendix B). The 
expansion for the flavor flelds given in Ea. (|A9|) leads to the following explicit form of the current operator: 

j^'(x,t) = // ^^e<^-^y-Y.\<i(*)^'<At)<eit)alAt)~v%At)r''vL^^^^^^ 

J J [2n)2 [2tt)2 ^,,1 
In Appendix O we prove that e,^t(0| J^(x, i)|0)e.;^ = 0. This result leads to (see Appendix A): 



1 1 



with 



^'(x,t) = A 



(2^)f 

Y^Ai) = sin2 0|C/k|xiX2 



7(k) 



- sm~ f |e """■•^"If/k 



1 



1 



where a ■ k 



and X I = 



We made use of the following relations: 

rt 



(18) 
(19) 

(20) 

(21) 
(22) 



a^J^Mj^S^^XuAty^'^^'S^k-p) , {/3l,,(t), a;t^(0)} = (<7.k)-yk.e(t)e-^"^'^*<S'(k-p)(23) 

The expression in Ea. (|19() contains the most general information about neutrino oscillations and can explicitly 
evaluated once the form of the wave-packet is specified. A similar expression can be easily obtained for the other 
quantity of interest, namely {i'e\J/^{'x.,t)\i'e)- 

For comparison with previous results and for better understanding the expression (|19|l . let us consider the limit 
situation in which the wave-packet becomes a plane wave, with definite momentum: /(k) = (2Tr)2S^{k — q). This 
obviously means that we loose information about localization and we need to integrate over the entire volume to get 
the flavor oscillations (in time). We indeed flnd (pe = J^): 



{iye\Qe{t)We) = / d^ K {ly^l Pe{x) We) = 1 - siu"^ (29) 



t I +|Vq|^sin^ 



2 ,-„2 / ^q,2 



(24) 



^ We work in the Heisenberg picture, the time evolution is thus borne by the operators. 
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which agrees with the result obtained in Refs.'^. Notice the presence of the non-standard oscillation term and of the 



momentum dependent amplitudes (the Bogoliubov coefficients satisfying \Uq 



We now consider the situation in which we have a spherically symmetric emission of the neutrinos from the source 
at (xo, to)- This allows us to limit our investigation to the radial flux, which can be identified without loss of generality 
with the z- component of the current. In this case, Ea. (|19|l takes it simplest form, since F'^ is diagonal in the chosen 
(chiral) representation. 

The matrix F'^ can be expanded as a linear combination of spinor outer products: 



10 
0-100 
0-10 
1 



= E I 



(25) 



where 























Vi = 
















Vo/ 





Vi = 



(26) 



This decomposition permits a remarkable rearrangement, since Ea. (|19|l then becomes: 



{V,\jl{-K,t)\v, 

where (see also Appendix ID|| 



E(-l)^+'*J(x,t) [W] *,(x,t) = E(-l)'''' lA' + ^ll 



(27) 



«'j(x,i) = ryt*(x,i) 



(28) 



Eq. H27|) gives a compact expression for the flavor current expectation value in the case of a spherically symmetric 
emission. It can indeed be integrated to give a useful expression for the electron neutrino flux at a given distance from 
the source. Although this calculation is of interest, because it contains information on the three-dimensional nature 
of the propagation, it will be not developed further here, and in the following Section we will consider the case of 
one-dimensional propagation. Note that in the present case (spherical symmetry), the expectation values of Jg and 
Jg at any point can be obtained from the above expectation value of J^. In the general case Ea. (|19|l . however, the 
expectation values of and have to be calculated independently and this can be done by use of the decomposition 
of F^ and F^ given in Appendix 

It remains to consider pe, which is the relevant measurable quantity for experiments without angular resolution. It 
is indeed easy to get its expectation value: 



(z^e|/Oe(x, t)\v^) 



which is related to the expectation value of as follows: 



(l/e|Je(x,i)ke 



(29) 



(30) 



This is a significant result in the sense that it directly contradicts the usual x ~ t approach: if the neutrino really 
travelled on a straight line then, for points for which the Cartesian z-direction coincides with the polar radial direction 
(as defined with respect to the source), one should have'^: (i/e| (a;)| i/g) = {vel p{x)\i'e)- 



^ the result would be the same up to a constant if one replaced the velocity of light by the (constant) average velocity of the neutrinos. 
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B. One— dimensional case 

Given the complexity of the three-dimensional analysis, we now limit the problem to one spatial dimension (as do 
most of other treatments in the literature) . 

The wave packet now takes the form /(k) = {2tt)^ S{ki) 6{k2) /(fcs), which corresponds to a packet with definite 
momentum in the x - and y - direction and a momentum distribution /(fca) along the z - direction. The uncertainty 
principle then dictates that the corresponding spatial distribution is that of a beam of infinite radius and uniform 
surface density moving in the z-direction (a neutrino sheet). 

In the one-dimensional case, we have (see Appendix IdIi : A2 — — and B2 = = for all r. In order to 
perform explicit evaluations, we choose in the following r = 1: this does not imply any loss of generality since the 
final result is independent of r. We also put Al + Bl = tp+ + tp'_^ and Al + Bl — t})- + 

Ea. l(?7|) then becomes: 



{jy,\J^{z,t)\i^e) ^ \^+{z,t) + i;'^{z,t)\ - \i;_{z,t) + ^pi{z,t) 



(31) 



with 



V'+(z,i) 
-0^(z,i) 
^_(z,i) 



A 



A 



dk 



Sin Ue ' 



Tfe(z) 



iA dk sin' d siniuk,2t)\Vk\ Tk{z)\Vk\ ~ T-kiz)\Uk\ 



dk 



cos2 6'e~"^'='i* ■ 



sin^ 9 e^"^"-^* 



T-fe(z) 



lA / dk sin^ 9 smiuk,2t)\Vk\ T-k{z)\Vk\ + Tk{z)\Uk\ 



(32) 
(33) 
(34) 
(35) 



where 



f{k) 



mi 



mi 



-ikz 



mi 



mi 



(36) 



The above relations are the main result of this work: starting with an arbitrary wave-packet f{k) we can now 
determine the expected one-dimensional flux of electron neutrinos at a given distance from the source as: 



$^^^^^(z)=/ dt{l^e\j!{z,t)\l^e 



(37) 



Note that the V'± terms are the corrections due to the use of the flavor Hilbert space in the definition of neutrino 
state. Ea. (|37l) is an exact expression for the space oscillation formula in ID and will be in the following evaluated 
numerically with the choice of a gaussian wave packet. In the following Section, however, we digress on the relativistic 
limit, which allows us to extract an (approximated) analytical form of Ea. l37|l and to compare it with the existing 
results in the literature. 



IV. APPROXIMATIONS AND NUMERICAL ANALYSIS 



In this Section, we explicitly evaluate Ea. H37|l by choosing a gaussian wave packet. We first consider the relativistic 
limit and show how the full result Eq. I|37|) reduces to the standard formula given in [25l l26l| . Then we proceed with 
numerical evaluations of the exact formula Ea. (|37l) . which allow us to test the non-relativistic regime. 

We choose f{k) to be a gaussian wave-packet with momentum spread ak- 



fik) = H exp 



(38) 



with H = {2TTal)~i. 
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A. Recovering the standard oscillation formula: sharp gaussian wave— packets in relativistic limit 



To recover previous results in the literature |25|,|26|, we consider the case of sharp gaussian wave-packet in relativistic 
limit, satisfying the following conditions : 

1. The mean momentum q of both mass species lies in the relativistic regime such that ujq^i ^ and \Vq\ — > 0. 
Hence ip'j^ and V'- can be ignored. 

2. The gaussian is sharply peaked, q 3> fft: the largest contributions to the integrals come from around the mean 
momentum q, i.e. Laplace's method |33| can be used. 

Using the above conditions, one has ( 1 ± ) " « 1 ± „"''' . Thus, to first order in 



ip+{z,t) !v AH dk 
Jo 



AH 



nil 



2uj. 



q.l Jo 



cos2 6'e~*("'=-i*-''^) + sin^eie"*^'^'''^*"'^^) 
dk cos'^ e 6-'^'^'"''+'''^ + sin2 6'e-*('^'=-^*+'=^) 





r (fc-g)n 


exp 


[ J 



expression: 



ip-{z,t) « AH 



mi 



dk 



cos^e'e-*^"'''!*-*'^) + sin^ee-'^'^'-'^*-''^) 





\ ik + q?] 


exp 


[ J 


we obtain for ip^ 




r (fc-g)21 


exp 


[ 4ai J 



(39) 



AH / dk 
Jo 



cos 





■ {k + q)^' 


exp 


[ 4a^ J 



(40) 



Again, the second term is exponentially suppressed and can be neglected. The first term of ip^ is also small for 
relativistic neutrinos and will be neglected in the following. We thus have: 

{l,e\j!iz,t)We) ^ \i^+iz,t)f (41) 



1. Momentum integration 



We now turn towards the evaluation of V'+ using Laplace's method. Let us first define the dispersion times 
rpdisp _ ^j^ .^2al, which set the timescale for the spatial dispersion of the wave-packets. 

Expanding the integrand about the mean momentum to second order in Sk — k — q, we obtain**: 



ip+{z,t) « AH cos^Oe 



-i{ujqAt—qz) 



exp 



(Z-Vg^tY' 



AH sin^^e-'^""'^*-*") exp 



4^1 (t) J.-, 



d(6k) exp 
d{Sk) exp 



-Si{t) [6k-i 



.{Z - Vq,it) 



-S2{t) [5k -i 



.(Z-Vg^t) 



252(0 



(42) 



4^2 (t) 

with Si{t) = cr^(l + it/Tf^^^). Also, is a spread in the configuration space, defined by (7j.ak = ^ and Vq_i = q/cuq^i. 

DO ^ pOO 

-q J — oo 

-Si' 



To a good approximation we can set ^ I^oo obtain 



oo 



d{6k) e 



1 + tt/T. 



(lisp 



Using the expansion a;^ ^ fa ujq^i + + Vq^iSk and making a change of variable k — q = 5k, we get 



1/-+ RiA/f cos^ee'f"^-'^'!'!*) / d(<5A:)e'*'=(^-"9,it) 



' exp 



-Sk^ 



1 . t 



9,1 



+AH sin2ee^(''^-"<J-2t) / e^'"=(^-"<J.2«) 



exp 



- 5k^ 



1 t 

T7 + i 

4(Tf. 2u), 



9,2 



which is brought in the given form by completing the square for the argument of the gaussian. 
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Following the usual treatment j^^, we now consider propagation for t ^ -j^disp^ Then Ii{t) becomes a constant: 
/ = 2y/Trak- For propagation times beyond J'^'^'p^ Laplace's method cannot be applied and one has to resort to the 
method of stationary phase 25, 33] . 

We are now in a position to calculate the interference term. We get, after some rearrangements. 



{iy,\j!{z,t)\iy,) ~ \A\^H^\I\^(^cos^9 e^p 



(Z ~ Vq.ltf 



2al 



sm V exp 



2al 



where the phase ^12 and jit) are given by 



■0sin2 0e^^^(*)-/(*) 



(43) 



(44) 



We note that ^(i) and /(f) are the usual terms obtained in the external wave packet model in QFT [251, |2g. We now 
proceed with time integration in order to finally get an expression involving distance only. 



2. Time average 

When calculating the interference term, spatial oscillations results from the cross-term. We show this explicitly by 
considering the first term in Ea. (|43|) : in the Laplace regime where [tlTf'^) « 0, if we assume that the experiment 
runs for a very long time. 



exp 



(Z - Vq^jtY 

2gI 



dt = 



2tt 



The first and second terms of Ea. H43l) are thus merely a collection of constants. 

Let us now turn to the oscillation term. The dominant contributions in the time integral come from exp[— /(i)], 
more precisely in the neighborhood of the maximum value of f{t), which occurs at 



Vq.l + Vq.2 



(45) 



Thus the integrand can be approximated as 



e-^m-m ^ exp 



i^max J max 



it 



dt 



2 dt^ 



(46) 



fmax ' f (irnax^ " 



«2 - <i) 



V2' 



-27r- 



z 



' coh 



(47) 
(48) 



to first order^. A coherence length = [L°^'^q) / {\/2'Kak) and an oscillation length = ATrq/Avn? with Am^ — 
•m\ — m\, having the same form as in the standard approach, are thus recovered. The other terms can be factorized 
as follows: 



4(T? 



^9,1+^9,2 



2ai 



2 \ y>q.l - ^q,2) 
"ll + Vq.2 



(49) 



^ We assume that f„,i ~ 1 



"12 « 2, vl. 



-t\m} jcp" and f5,ifg,2 ~ 1- 
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The first term integrates to a constant: yph^a^. As for the second term, it becomes '^(uj^ x — lOq.-i)^ ~ 27r^ (7%^)^ 
which gives the locahzation termj25|. Hence we can finally write down the (normalized) highly-relativistic space 
oscillation formula (compare to Ref . |2q ) : 



$,^^,^(z) ~ l-isin2(20)|l-cos(2^-^) exp 



" coh 



2'K- 



(50) 



Thus we have shown that the exact formula Eq.(|37|) reproduces the usual oscillation formula in the relativistic limit. 

B. Numerical Calculations 

In this Section, we perform numerical evaluations of the exact oscillation formula Ea. H37|l in order to compare it 
with the approximate result Ea. (|50|l for sample values of the parameters. We plot the two expressions in the maximal 
mixing case for a given au and for different values of q and of the masses: in the relativistic case (Fig.l) we observe a 
perfect agreement of the two formulas, as expected. 

We then explore the non-relativistic region and observe that relevant deviations from the standard formula do 
indeed appear, both in the oscillation amplitude and in the phase. This is already evident in Figs. 2, 3 where we use 
q = 100 with rrij = 1, = 3 in the first case and — 1, = 10 in the second case. The fact that the two plots 
are simply scaled with respect to each other, indicates that for these values of parameters the usual relation coherence 
length vs. oscillation length: L'^°^ = {L'"''^q) / {y^nak) is still valid with = A-Kq/Am'^. In Fig.4 we use q = 50 
with — 1, 771,2 = 3 and as expected we observe a larger deviation from the usual formula. 

Note that the value of the Bogoliubov coefficient \Vk\ is very small in all the considered cases®, so the observed 
corrections originate from the "standard" terms ip± Eos. (|32|l . (13411 in the oscillation formula Ea. (|37|) rather than from 
the flavor vacuum contribution, which thus turns out to be very difhcult to detect when considering (time-) averaged 
oscillation formulas like in the present case. 




Figure 1: Plot of the QFT flavor flux (solid line) against the standard oscillation formula (dashed line) 

for e = 7r/4, CTfe = 10, = 1, = 3, q = 1000. 



The values of \Vq\ are; 1 x 10"^ (Fig.l), 1 x IQ-^ (Fig.2), 4 x lO'^ (Fig.3), 2 x IQ-^ (Fig.4). 



■] , , , . . , 

200 400 600 oUU lUUO 1200 

Figure 2: Plot of the QFT flavor flux (thick hne) against the standard osciflation formula (light line) 
for 6 = 7r/4, ak = 10, = 1, = 3, q = 100. 



1 
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Figure 3: Plot of the QFT flavor flux (thick line) against the standard oscillation formula (light line) 
for 9 = 7r/4, ak = 10, = 1, = 10, q = 100. 
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Figure 4: Plot of the QFT flavor flux (thick line) against the standard oscillation formula (light line) 

for 6 — it/4:, at = 10, = 1, = 3, q — 50. 

A more complete discussion of these corrections and of their possible phenomenological relevance will be given 
elsewhere ,34j . 



V. CONCLUSIONS 



In this paper, we have for the first time derived a space-time dependent oscillation formula directly from the 
relativistic currents for flavor fields. These currents have been recently studied in Ref.^3]) a result which served as 
basis for the present analysis, together with previous results on the Hilbert space for mixed fields'^, 

We first presented a general expression for the electron neutrino fiux in three dimensions and then specialized to 
the case with spherical symmetry, for which we were able to find a more explicit expression. In order to perform 
further analysis and numerical evaluations, we then considered the one-dimensional case with gaussian wave-packets, 
which is also the one most frequently treated in literature. 

Our formulation presents several advantages with respect to existing treatments of neutrino oscillations in Quantum 
Field Theory: it is a very straightforward approach which is easy to relate to practical experimental situations; it 
takes into account the non-trivial nature of the flavor vacuum and flavor states are thus consistently defined. It 
takes explicitly into account the full spin structure of neutrino states and does not resort to relativistic limit and/or 
assumption of nearly degenerate masses for the energy eigenstates. 

We have shown how, in different limits, our formula reproduce existing results. Thus, in the case of relativistic 
neutrinos with nearly degenerate masses, we recover analytically the standard space dependent expression for neutrino 
oscillations |25y2Q], which is thus once again confirmed from an independent approach. Also, previous results on the 
flavor chargesjoj have been recovered, exhibiting the non-standard oscillation terms. 

The numerical analysis shows that in the non-relativistic regime, our formula predicts significant deviations from 
the standard oscillation formula 25, 26], which on the other hand cannot be expected to be valid in that region. An 
analysis of the phenomenological implications of the results obtained in this paper will be presented elsewhere js^. 
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APPENDIX A: FLAVOR HILBERT SPACE 



The free fields J^i(a;) and 1^2(2^) are written as {t = xq) 



vAx) = >^ , „ 



i = 1,2 



(Al) 



where ul^^{t) = e-^'^'^-^'u^; and vl,^(t) = e*'^'=-*<,, with = ^/W+^^■ The a^,^ and the /3^^,, i,r = 1,2 are 
the annihilation operators for the vacuum state |0)i^2 = |0)i ^ \^)2'- aki|0)i2 = /3ki|0)i2 = 0. The anticommutation 
relations are: 



All other anticommutators are zero. The orthonormality and completeness relations are: 

r 

where In is the n x n unit matrix. Eq.(|2|l can be recast in the form: 

i^,{x) = G-\t)v,{x)Ge{t), 

Gg(t) — exp 9 j d'^-x.(^\{x)v2{x) — i'\{x)i>i{x)^^ 

with (cr, i) ~ (e, 1), (/x, 2). The generator Gg{t) does not leave invariant the vacuum |0)i^2: 



2 • 



(A2) 
(A3) 

(A4) 

(A5) 
(A6) 

(A7) 



We will refer to |0(t))e,p as the flavor vacuum: it is orthogonal to |0)i^2 in the infinite volume limit,S|. We define the 
fiavor annihilators, relative to the fields Ve{x) and Vi_i{x) as^ 



alAt)^G,\t)al,(t)Ge{t) ; pX-^{t) ^ G,\t)(3X^{t)Ge{t) 



with (cr, i) — (e, 1), (/x, 2). The flavor fields can be expanded as: 

.=1,2^ (27r)^ 



i^<T(a;) = E 



(A8) 



(A9) 



with (cr, i) = (e, 1), (/i, 2). The flavor annihilation operators are defined as operators 





— COS t 


' «Li - 


- sin 6* 


s 




(AlO) 




= cos ^ 


'«L2 - 


- sin 6 


EK2W<iW«k,i + <2W^-k,iW/3l,i 

s 




(All) 




= COS^ 


' /31k,i 


+ sin 


E [A,2W«:k,lW /31k,2 + <2W«!lk,lW <2 

s 


(A12) 




= COS ^ 


' /3:k,2 


— sin 


^ E [A,iW«1k,2W /3ik,i + <iW«1k,2W <i 


(A13) 



The annihilation of the flavor vacuum at each time is expressed as: e(i)|0(t))e,fi = Gg ^(t) ^ |0)i,2 = 0. 



In the reference frame where k is colUnear with k = (0, 0, 1), the spins decouple and we have: 

aljt) = cos9al, + sinO {u^{t) a^^ + eiV^it) p%) 
Q^k.^C*) = cos6'Q;k2 - sinfi" (?7k(i) - 4 ^k(i) /3!!.k_i^ 
/3-k,eW = cosepLu,! + sine (uUt)f3Lu,2 - elVi,{t)al 
/3-k,^(i) = cose/3:k,2 - sin^ (c/kW Ck.i + 4 Vic(t) <i) 
where = (^_2^y+k k+i g^j^^j Uu{t), Vk(i) are BogoHubov coefficients given by: 



satisfying |C/k|2 + |ykP = l. 



+ mi \ ^ / ^^fe,2 + \ / Ik| |k| 



2wfe,i / V 2a;fe,2 / V(wfc,2 + m2) (wfe,i+mi) 



APPENDIX B: CHIRAL REPRESENTATION AND USEFUL RELATIONS 



ujk,i+mi 
•k 



, r = l,2 



It 1 2+ 2 o / , 
^-k,l^k,2 = -V-kAUk,2 = 2XlX2 — — + 



V-k 1^1,2 = (^'-fe i^m) =2x1X2 ( ^ + 

and u^\ui. .j = for r 7^ s. We define: 



Wfc,i+mi Wfe,2+m2, 

—) 

2 + m2 / 



ro = i4 , r^ = ( ) , » = 1,2,3. 



Notice that and can be decomposed as foUows: 



with 




Ai = n > ^1 
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APPENDIX C: VACUUM EXPECTATION VALUE OF THE CURRENT 



We show here that e,/i(0| J^(x, t)|0)e.p — 0. Let us consider for example the foUowmg quantity: 



e,/j \ 



«p!eW<e(i) |0)e,^ = 1.2 (0 1 (0) t)a^,, (i)G,- ^0) 



11.2 



(CI) 



Now let us define aj^ ^{t) = G6/(0) aj^ ^{t) Gg ^(0) and h.c. Of course, ^(0) = a\^i- It is easy to realize that (see 
Eqs.Eni)): 

{5LeW,5^!e(0} - for k^p, Vt,i',r,s. (C2) 

This in turn implies that the quantity l|(yl() vanishes for k 7^ p. Since this is valid for any two flavor operators, we 
only need to calculate the various terms for the (flavor) vev of the current for equal momenta. We thus write: 



(0|J,^(x,i)|0)e,^ = / d=^ke,^(0|J,^(k,i)|0)e,^ - 0. 



(C3) 



The flavor vacuum is invariant under rotations, so we can choose without loss of generality the most suitable spatial 
reference frame for each of the terms in the momentum integration. In our case, this is the one for which k is collinear 
to (0, 0, 1). In this frame, we obtain the relations: 

l/3:leW/3ik,eW|0)e,M 



e.M(0|<eW«k,eW|0) 



e,/.i 



cos2 e sin^ ( ^^l±^t ) + sin^ 0|C/kp sin^ {cv^t) 



(C4) 



V 2 

e.M(0|/3:k,eW«k,eW|0)e,p = ( e.p (0 1 a;,^ t)/?:^^, (<) 1 0) 

= S''' el.sm'^ e^"^^* |C/k| |Vk| (l + e-"^"^^^ - 2e-2H'^i+'^2)t^ cos^ 6* +2isin{uj2t) (e'^'^^'lt/kP + e*""*| VkP) ain\m:>) 

Because of relation (jC4|) and the orthonormality relations for the spinors, we easily realize that e.^i(0| (k, t)\0)e,^ — 0. 
As for e,^i(0| Jg (k, i)|0)e,;j, we need to consider: 



„'-t 



(C6) 
(C7) 



As a consequence of relations (|C4|I . (|C5|I and ljC6|l . ljC7|l . the expectation value e,p(0| Jg (k, i)|0)e,^i is an odd function 
of k and therefore its integral vanishes. 

Note that e.^^{0\J^''^{i^,t)\0)e.^ vanish identically in the chosen reference frame, because of: 



(C8) 



APPENDIX D: EXPLICIT FORMS 



^^(x,<) = A 



(2^ 



J (27r)2 L 



<5^"Xk.e(i) 

(a • k)^'^Yk.e(i) 



(Dl) 
(D2) 



/d^\c f d^V 

— e*- /(k) XI (1 + fca/f^i) J^k,e(i) ; = ^ / -— T e*-'' /(k) xi (^3 - Ikp/^^i) Y^^A^) 

(27r)2 J (27r)2 

— e^'^-- /(k) XI (fc+/f^i) Xk,e(t) ; =-Bl=A — e''^-'' /(k) xi k+Y^At) 

(2n)2 J {2n)2 



- -Al = A 



— e'''-- /(k) XI (1 - fca/f^i) Xk,e(i) ; Bl = -A — e^'^ '^ /(k) xi (|k| V^^i + ^^3) ik,e(t5D3) 
(27r)2 J (27r)2 
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where fii = w^.i + mi and = ki + ik2 , fc_ = fci — ik2. 

For k coUinear with k = (0, 0, 1) or in the one dimensional case, we have 

{/3lVeW, <e(0)} = '5^''ek|f/k||Vi| sin^ 6 [e'^^-^' - e"'-'-^*] e-'"^'^* . (D4) 
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